A spatially shifted beam approach to subwavelength focusing 
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Although negative-refractive- index metamaterials have successfully achieved subwavelength focus- 
ing, image resolution is limited by the presence of losses. In this Letter, a metal transmission screen 
with subwavelength spaced slots is proposed that focuses the near-field beyond the diffraction limit 
and furthermore, is easily scaled from microwave frequencies to the optical regime. An analytical 
model based on the superposition of shifted beam patterns is developed that agrees very well with 
full-wave simulations and is corroborated by experimental results at microwave frequencies. 
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Subwavelength focusing of electromagnetic radiation 
has been the subject of a great deal of recent research 
in disciplines ranging from the microwave to the optical 
regimes. When electric and magnetic fields arc produced 
or scattered, the decay of evanescent field components 
away from their source or scatterer leads to a loss of spa- 
tial resolution below the order of one wavelength. Many 
efforts have been made to overcome this effect, known 
as the diffraction limit, with various methods achieving 
a measure of success. Early efforts based on the use of 
small apertures were very successful 1], leading to the 
field of near- field scanning optical microscopy (NSOM), 
but required extremely precise probe placements that in- 
creased the cost and complexity of the techniques. In 
2000, J. B. Pendry proposed an alternative strategy with 
a "perfect lens" that could focus electromagnetic waves 
to a point through the amplification of evanescent field 
components This was successfully demonstrated ex- 
perimentally in planar form Q and in the quasistatic 
limit [1, 0| but both cases suffered greatly from losses. 
Other methods that have seen some recent developments 
are based on field synthesis through the plane-wave ex- 
pansion of aperture fields and include direct field inte- 
gration leading to near-field plates and transmission 
screens inspired by holography fi\ . Here we show a differ- 
ent conceptual approach to synthesizing such transmis- 
sion screens that is both simple and intuitive, incorporat- 
ing the superposition of properly weighted dipole slot an- 
tenna field patterns. The key enabling observation is that 
the beam patterns formed by displaced slot elements are 
spatially shifted in the near-field, allowing them to form 
a set of basis functions for a moment expansion synthesis 
of image-plane field distributions. This technique leads 
directly to physically realizable implementations and ex- 
perimental validation that is not sensitive to loss and can 
be readily scaled up to optical frequencies. 

Field patterns with subwavelength features contain 
spatial Fourier components with transverse propagation 
vectors kxz larger than k = uj^yeJI. To satisfy Maxwell's 
equations, these components decay evanescently along y, 
and are therefore restricted to the near-field. Negative- 



FIG. 1: A qualitative illustration of how an array of dipole 
slots produces shifted beams in the near-field but identical 
magnitude distributions in the far-field. 



refractive-index slabs can be used to recover the near- 
field through evanescent field amplification but the im- 
age resolution tends to be very sensitive to losses 0, 3 . 
Without such lenses, the image plane must remain close 
to the source since the field profile broadens quickly to 
the diffraction limit. In this Letter, closely spaced dipole 
slots less than half a wavelength apart produce shifted 
beams that are weighted with out-of-phase currents to ef- 
fectively increase the range in which the evanescent fields 
can be used to perform subwavelength focusing. 

The electric field produced by an infinitesimally thin z- 
directed dipole slot antenna located at the origin is given 
in Eq. ([1]) for a range of fields restricted to the xy-plane 
Q. The slot is fed by a magnetic source current /„, with 
the amplitude of the electric field related to the distance 
p from the center of the dipole and the phase related to 
the distance R from the tip of the dipole. In accordance 
with the restrictions imposed by the experimental setup, 
only the x-component of the electric field is considered 
in this Letter, but any other component of the electric or 
magnetic field can be focused using the same principles 
outlined for E^.. 
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It is well known in optics and electromagnetics that 
the far-field distribution of a radiating system is simply 
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the Fourier transform of the source field distribution Q . 
This implies that if a slot were moved from the origin to a 
new position x = d along the x-axis, the far-field pattern 
would experience only a phase shift, as shown in Eq. (|2ap . 
The separation of the far-field into the product of an el- 
ement factor given by Eq. ([T|) and an array factor is not 
valid in the near-field, however. Indeed, the situation 
is distinctly different, giving us a new handle to manip- 
ulate the near-field for effects such as subwavelength fo- 
cusing. The near-field distribution of an antenna element 
translated along the x-axis is the spatially shifted beam 
described by Eq. (j2b[) . The presence of strong evanes- 
cent waves produces fields that decay rapidly away from 
the slot element, leading to a localization of fields shifted 
away from the array origin. This is critical in providing 
sufficient amplitude variation to design a subwavelength 
focused beam. As the distance from the source increases, 
the shifted beams spread until they can no longer be dis- 
tinguished from one another strictly on the basis of am- 
plitude and the property of spatially shifted fields is lost, 
as demonstrated in Fig. [T] 



E:,,sMfted = E,{x,v,z)e-^^''-°^'^ (far field) (2a) 
Ex,shifted = Ex{x~d,y,z) (near field) (2b) 



For a source to focal plane separation of y = 0.15A, 
the individual dipole fields from a three-dipole array are 
plotted in Fig. [2^. By adjusting the complex weights 
on each of these element field profiles (Fig. [2]d), the two 
satellite slots can be designed to destructively interfere 
with the central profile and produce a null in the total 
field close to the peak of the main lobe (Fig. [2t). As 
the distance between the nulls decreases, the main lobe 
beam width becomes smaller and the beam is more fo- 
cused. A tradeoff is involved, however, since the element 
beam profiles are of finite width themselves, leading to 
increased side-lobe levels and higher main lobe attenua- 
tion relative to the constituent element fields. Increasing 
the number of elements can serve to lower the side lobes 
but a limit is reached as to the minimum focused beam 
width for a given inter-element spacing d. 

The analytical formulation for the beam produced by 
a single slot can be extended to a linear array with an 
arbitrary number of slot elements. The total x-directed 
electric field for an array of N elements is given by Eq. ([3]) 
where Wnim is the n*'' weighted magnetic source current 
and X — nd\s the position along the x-axis of the n*'' slot. 
In the vector representation, w is a iV-element column 
vector containing the complex weights w„ assigned to 
each slot and Ex is a A'^-element row vector containing 
the field from each element as a function of x. 
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FIG. 2; (color online) Subwavelength focusing from the su- 
perposition of shifted beams in the near-field. The transverse 
electric field at y = 0.I5A is plotted at several conceptual 
stages, (a) Shifted beams from slots at a; = —0.1 A (red), 
X = Q (blue), and a: = 0.1 A (green), (b) The shifted beams 
multiplied by their weighting factors, (c) The weighted shifted 
beams superposed to produce a subwavelength focused beam 
(solid) and the beam from a single slot (dashed). 
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The weights Wn are chosen to scale the individual slot 
electric field components E^^n to produce the desired 
electric field profile at the focal plane. In order to de- 
termine the optimal weights Wopt necessary to approx- 
imate a desired electric field profile f{x), the individual 
slot field components can be regarded as basis functions 
for a method of moments field expansion. By multiply- 
ing both sides of f{x) ~ Ex • Wopt by E* and integrat- 
ing over X we are left with a system of N independent 
equations that can be solved for the N weights Wn by 
reducing the problem to that of inverting an N x N ma- 
trix. Although these basis functions are linearly indepen- 
dent, they are not orthogonal so the off-diagonal terms 
of the matrix do not vanish and the analytic formula- 
tion becomes unwieldy. Instead, f{x) and E^^nix) can 
be written as column vectors b = /[x] and A = -E^.nia^] 
by discretizing them over x, causing the method of mo- 
ments to be reduced to the method of least squares 
with an expression for w given in Eq. ^ . 
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The desired focusing effect is demonstrated using 
a Gaussian pulse with a full-width half-maximum 
(FWHM) beam width of 0.12A as the target field profile 
for an array of three slots separated hy d = O.IA. The op- 
timal weights were calculated using Eq. ([4]) to be tuo = 1 
and ■w±i = 0.446Z174° for a focal plane at y 0.15A 
(Fig. [2fc). 

In order to show the increased focusing power of an 
array with more elements, the focal plane was moved 
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FIG. 3: (color online) Simulated focused beams for slot ar- 
rays with 3 and 9 elements optimized for a focal plane at 
y — 0.25A. The thin red line was calculated analytically from 
Eq. ([si while the thick blue line was extracted from a full- 
wave simulation. The dotted line indicates the beam formed 
from a single slot, (a) Three-element array with weights 
Wo = 1 and w±i — 0.497/172°. (b) Nine-element array 
with weights wo = 1, it)±i = 0.764/179°, w±2 = 0.373/-2°, 
w±3 = 0.136/175°, and w±4 = 0.032/-6°. 



to y = 0.25A where the beam width of a single slot is 
0.5A. The optimal weights were calculated for an = 3 
and an = 9 array with d = O.IA. The structure was 
simulated using Ansoft's HFSS full-wave finite-element 
software with the half- wavelength dipole slots fed by ideal 
x-directed voltage sources. The resulting field profiles 
matched very well with those calculated using Eq. ([3]), 
converging to the analytical solution as the width of the 
simulated dipoles was reduced. The three-element array 
was able to produce a beam width of 0.250A (Fig. [3^) 
while the nine-element array focused the beam further to 
a beam width of 0.144A (Fig. Oj). 

In a numerical implementation, Eq. (j4]) must be calcu- 
lated using a truncated range of x values, typically chosen 
such that the target field is very small at each endpoint. 
When the weights are optimized over a small portion of 
the target function, however, such as by truncating the 
range of x to those values falling within the target func- 
tion's FWHM beam width, an arbitrarily narrow central 
beam can be formed for any iV > 3 and any d. However, 
this comes at the expense of extremely large side lobe 
amplitudes and high weight sensitivity, an observation 
reported in super-directive and super-oscillatory studies 
in the literature that can render ambitious applications 
of these procedures impractical [ll|, • If the weights 
are instead optimized over several wavelengths in x, the 
side lobes are reduced to reasonable levels but the narrow 
beam width is no longer maintained. Increasing the num- 
ber of elements reduces the side lobes further and allows 
a narrow beam width at more distant focal planes. 

In order to realize this structure for a practical exper- 
iment, the ideal voltage sources at the centers of every 
slot were removed and plane wave incidence was consid- 
ered. As demonstrated through full-wave simulations in 
0, the weights were implemented by changing the di- 
mensions of each slot, with a narrower slot producing a 
higher amplitude electric field and a wider slot producing 
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FIG. 4: (color online) Simulated electric field magnitudes over 
the xy-plane for a single-slot (a) and a triple-slot (b) array 
under plane wave excitation at 10 GHz. The fields are nor- 
malized along the i-direction with respect to the field value 
at a; = to highlight the profile of the beam. The thin black 
contour trace indicates the FWHM beam width. 
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FIG. 5: Full-wave simulation results for the magnitude of the 
electric field along the main axis of the focused beam as a func- 
tion of distance from the screen. The dashed line indicates 
the field strength of the incident electric field while the solid 
line indicates the field strength of the focused beam. At the 
focal plane {y = 0.15A), the focused electric field strength is 
larger than the incident field by a factor of 1.54. Measurement 
points from experiments conducted at j/ = O.IA, y = 0.15A, 
and y — 0.2A are included, showing a general agreement with 
simulations. The magnitude of the measured field is normal- 
ized with respect to the background electric field (measured 
with the screen removed) and has a value of 1.33 at y = 0.15A. 



a lower amplitude electric field. Implementing the near 
180° phase shift between adjacent slots was more chal- 
lenging and necessitated increasing and decreasing the 
lengths of adjacent slots to produce inductive and capac- 
itive slots that could resonate together. The final slot di- 
mensions were optimized iteratively using full-wave sim- 
ulations until the desired total field pattern was achieved. 
Fig. 2^ and Fig. |3|3 present a cross-section of the fields 
produced by the slots for a single-slot and a triple-slot 
arrangement, respectively. 

Although not visible in the normalized data of Fig. 
m despite the presence of rapidly decaying evanescent 
fields, the focused electric field magnitude at the focal 
plane is 1.5 times the magnitude of the incident field due 
to slot resonances (see Fig. [5]). The fields also show a 
similarity to those observed in negative-refractive-index 
slab focusing (see Fig. 3 and Fig. 4 from [3]). The 
presence of evanescent field components in any type of 
near-field focusing implies that rather than producing a 
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FIG. 6: (color online) (a) Experimental measurements of the 
transverse electric field magnitude (solid lines) for single-slot 
and triple-slot arrays at 10 GHz compared with full-wave sim- 
ulations (dotted). The slight asymmetry in the measured field 
profile resulted from the asymmetrical coupling between the 
slots and the probe antenna (see inset), (b) Measured field 
magnitudes at three frequencies for the single slot (dotted) 
and triple slot (solid) screens. 
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FIG. 7: Simulation results from a full- wave simulation of 
the three-slot focusing screen adapted for optical frequencies 
(A = 830 nm). The transverse electric field is plotted at the 
focal plane O.IA away from the screen. The metal screen was 
implemented using a lossy Drude model for silver with a fi- 
nite thickness of 40 nm. The slot lengths were shortened to 
200 nm and 130 nm to compensate for the increased capaci- 
tance of slots at plasmonic frequencies. The focused FWIfM 
beam width was reduced to 63% of the single slot beam width. 



three-dimensional spot, the fields are confined only in the 
transverse direction One major advantage of the 

focusing screens over metamaterial lenses, however, is the 
fact that Ohmic losses occm only at the screen and not 
throughout a region of wave propagation. The structure 
is therefore insensitive to loss and although much of the 
incident power is reflected at the screen, this in no way af- 
fects the focusing behavior. The simplicity of the slotted 
screen is also a distinct advantage over bulk metamate- 
rials, allowing it to be easily scaled in frequency. 

The dipole slots were created by cutting a three-dipole 
pattern in a 6 mil stainless steel sheet. The slots were 
designed to operate at 10 GHz with a central slot dimen- 
sion of 13.2 mm x 1.2 mm and a satellite slot dimension 
of 17 mm X 0.6 mm spaced 3 mm away from the central 
slot. An X-band horn antenna with a Rexolite bicon- 
vex lens was used to generate a collimated beam. The 
slot apertures were placed at the Gaussian beam waist, 
forming a very good approximation to a plane-wave. An 
Agilent E8364B network analyzer collected the field data 
over 8-12 GHz with one port connected to the horn an- 
tenna and the other port connected to a co-axial cable 
probe whose exposed inner conductor was bent ninety 
degrees to form a miniature dipole receiver antenna (see 
Fig. [6^ inset). A custom built XYZ-stage from Newmark 
Systems positioned the probe over the screen with high 
precision to collect the near-field data on a plane 4.5 mm 
above the slotted sheet. 

To verify the experimental setup and calibrate the 
probe, the fields produced by the horn and lens setup 
were measured with the screen removed and found to be 
collimated at the screen plane. The transverse electric 
field was measured over the slots and the focused beam 
width was found to be 5.2 mm, which was within 6% of 
the simulated value of 4.9 mm. The two satellite slots 
were then covered with 3M copper tape and the mea- 
surements were retaken, resulting in a single-slot beam 



width of 9.7 mm, within 1% of the simulated value of 9.6 
mm. The measured beam width was reduced by a factor 
of 1.87 through the use of the two satellite slots (see Fig. 
|6^). In Fig. [6)3, two other frequencies close to the design 
frequency are plotted to show the frequency dependence 
of the focusing screen. 

In summary, a shifted-beam approach to near-field fo- 
cusing using slot antennas has been demonstrated in this 
Letter. Using an array of dipole slots weighted out-of- 
phase with respect to one another, a set of shifted beams 
is formed that yields a subwavelength focused beam 
through superposition. A moment-expansion method 
was developed to find the optimal weights given the 
number of elements, the element spacing, and the tar- 
get field profile. Out-of-phase weights were implemented 
for plane-wave excitation using resonating inductive and 
capacitive slots. This led to experimental verification of 
the shifted-beam method for an iV = 3 array, with tighter 
focusing observed in simulations using larger arrays. 

The inherent scalability of the slot apertures provides 
the opportunity to implement the focusing screen at fre- 
quencies ranging from microwaves to Terahertz and all 
the way up to the optical domain for corresponding ap- 
plications such as high-resolution imaging and sensing. 
The dependence of the dipole slot behavior on plasma 
oscillations present in metals at optical frequencies can 
be readily accounted for by introducing an effective slot 
length that restores the scaled resonance frequency 
This was successfully demonstrated in simulations at 
A — 830 nm with slots cut into a 40 nm silver slab, 
as shown in Fig [71 The increased subwavelength focus- 
ing range is very attractive for near-field imaging appli- 
cations which are often constrained to working within 
an extremely small fraction of a wavelength away from 
the object to be imaged. The structure demonstrated in 
this Letter can be used to resolve multiple objects in a 
variety of two-dimensional configurations with subwave- 
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length precision by mechanically scanning the structure 
using orthogonal scanning patterns and detecting the 
back-scattered fields. For complete imaging, the struc- 
ture can be extended fully to two dimensions by adding 
a second set of satellite slots offset slightly along x and 
shifted along z. 
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